Spin-triplet Supercurrent through Inhomogeneous Ferromagnetic Trilayers 



Mohammad Alidoust 1 and Jacob Linder 1 

'Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

(Dated: February 19, 2011) 

Motivated by a recent experiment [J. W. A. Robinson, J. D. S. Witt and M. G. Blamire, Science, 329, 5987 
(2010)], we here study the possibility of establishing a long-range spin-triplet supercurrent through an inhomo- 
geneous ferromagnetic region consisting of a Ho|Co|Ho trilayer sandwiched between two conventional s-wave 
superconductors. We utilize a full numerical solution in the diffusive regime of transport and study the behavior 
of the supercurrent for various experimentally relevant configurations of the ferromagnetic trilayer. We obtain 
qualitatively very good agreement with experimental data regarding the behavior of the supercurrent as a func- 
tion of the width of the Co-layer, Lq - Moreover, we find a synthesis of 0-7T oscillations with superimposed 
rapid oscillations when varying the width of the Ho-layer which pertain specifically to the spiral magnetization 
texture in Ho. We are not able to reproduce the anomalous peaks in the supercurrent observed experimentally 
in this regime, but note that the results obtained are quite sensitive to the exact magnetization profile in the Ho- 
layers, which could be the reason for the discrepancy between our model and the experimental reported data for 
this particular aspect. We also investigate the supercurrent in a system where the intrinsically inhomogeneous 
Ho ferromagnets are replaced with domain- wall ferromagnets, and find similar behavior as in the Ho | Co | Ho 
case. Furthermore, we propose a novel type of magnetic Josephson junction including only a domain-wall fer- 
romagnet and a homogeneous ferromagnetic layer, which in addition to simplicity regarding the magnetization 
profile also offers a tunable long-range spin-triplet supercurrent. Finally, we discuss some experimental aspects 
of our findings. 

PACS numbers: 74.45,+c, 74.50.+r, 75.70.Cn, 74.20.Rp, 74.78.Na 



1. INTRODUCTION 



Due to the rich physics from a fundamental view- 
point and the possibility of birthing practical applications, 
configurations containing superconductors and ferromag- 
nets have attracted much attention theoretically 4-8 ' 11,12 and 
experimentally 1-3 ' 9 ' 10 ' 13-15 over the last years (see for instance 
Ref. 2 for a comprehensive review and reference-list). One 
of the most intriguing phenomena in the context of this in- 
terplay is the generation of a long-ranged spin-triplet super- 
current flowing through a Josephson junction with magnetic 
elements. The main criterion for generation of such a cur- 
rent is that some form of magnetic inhomogeneity must be 
present in the junction, as first shown by Bergeret et al. 20 ~ 22 
and Volkov et al. 23-25 . In the diffusive limit of transport, being 
often the experimentally most relevant regime, such a non- 
uniform magnetization can induce exotic long-range super- 
conducting correlations which are odd under time-reversal: 
the so-called odd-frequency superconducting state. The re- 
sulting triplet-supercurrent decays over same length scale as in 
a superconductor|normal metal | superconductor (S|N|S) junc- 
tion, but now additionally activates the spin-degree of free- 
dom in the transport of Cooper pairs. A long-range triplet- 
supercurrent was predicted to occur in a setup consisting 
of three non-collinear homogeneous magnetic layers 15 , but 
disappears in the scenario of only two homogeneous non- 
collinear magnetic layers, as discussed in Refs. 26,35,38,39. 
Over the last couple of years, induction of triplet-correlations 
in layered heterostmctures with ferromagnets (F) and super- 
conductors (S) in the clean limit has been studied by using 
different formalisms and configurations in Refs. AQ-A2. On 
the experimental side, Keizer et al. 44 observed a long-range 
supercurrent through half-metallic Cr02, whereas very recent 



work also reports observation of a long-range supercurrent 
through an inhomogeneous magnetic layer 6 ' 27 . In particular, 
Robinson et al? 5 investigated the appearance of a spin-triplet 
supercurrent flowing through a magnetic Ho|Co|Ho trilayer. 
Due to the intrinsic magnetic inhomogeneity in Ho, featuring 
a spiral magnetization texture, it was found that a strong spin- 
triplet supercurrent was established through the trilayer con- 
necting two s-wave superconducting leads. In this paper, mo- 
tivated by the very recent experimental in Ref. 35, we utilize 
a full numerical solution of the quasiclassical Green's func- 
tion in the diffusive regime and study theoretically spin-triplet 
condensation in the critical charge current flowing through a 
Ho|Co|Ho magnetic trilayer. Due to our numerical approach, 
we have access to the full-proximity effect regime and com- 
plicated magnetization textures in the trilayer. This allows us 
to also study the influence of domain-walls in the ferromagnet 
on the behavior of the long-range supercurrent. 

The main results in Ref. 35 due to Robinson et al. were (i) 
a slow decay of the supercurrent as a function of the Co-layer 
thickness and (ii) anomalous peaks arising in the characteristic 
voltage of the junction as a function of the Ho-layer thickness. 
Using the computational machinery described above, we ob- 
tain qualitatively very good agreement with the experimental 
data pertaining to (i). However, we are not able to reproduce 
the anomalous peaks observed for the supercurrent pertaining 
to (ii). Instead, we find a synthesis of Q-tt oscillations with su- 
perimposed rapid oscillations which pertain specifically to the 
spiral magnetization texture in Ho. However, we also show 
how the exact behavior of the supercurrent vs. the width of 
the inhomogeneous magnetic layer is rather sensitive to the 
exact magnetization pattern. This suggests that the trilayer 
magnetization texture realized in the experiment by Robin- 
son et al? 5 might differ somewhat from our model. Moti- 
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vated by the above mentioned reason, we also investigate how 
the critical current behaves when we replace the Ho layers 
in the trilayer junctions with domain-walls. Finally, we pro- 
pose a novel type of inhomogeneous ferromagnetic Joseph- 
son junction with a simpler magnetization profile compared to 
previous proposals, including only a domain-wall and homo- 
geneous ferromagnet, and demonstrate the possibility to tune 
the long-ranged spin-triplet supercurrent flowing through the 
junction. 

This work is organized as follow: In Sec. , we present 
the main ingredients of the theory which is used throughout 
the paper, i.e. a quasiclassical Green's function method in 
the diffusive limit studied by means of the Usadel equation 
and supplemented with proper boundary conditions. In Sec. 

, we investigate the triplet-supercurrent flowing through an 
Ho | Co | Ho ferromagnetic trilayer as a function of both the Ho- 
and Co-layer thickness including two different magnetization 
textures in the Ho layers. In Sec. , we investigate how the be- 
havior of the supercurrent is altered when the Ho regions are 
replaced with domain-wall ferromagnets, which also feature 
an intrinsical inhomogeneous magnetization texture. We also 
propose a novel type of ferromagnetic Josephson junctions to 
investigate the possibility of tuning the long-range contribu- 
tion to the supercurrent via an external field. Finally, we sum- 
marize and give concluding remarks in Sec. 5. 

2. THEORETICAL APPROACH 

For studying various characteristics of different media, the 
Green's function method is a fundamental approach utilized 
in many areas of condensed-matter physics 29 ' 30 . In the con- 
text of non-equilbrium transport through different media one 
should start from Dyson's equation of motion and calculate 
the Keldysh Green's function 31 . In equilibrium situations, 
the Keldysh block of Green's function can be obtained from 
the Advanced and Retarded blocks of the Green's function. 
Inside superconducting regions, Dyson's equation of motion 
equation transforms to Gorkov's equations which in turn can 
be reduced to Eilenberger's equation within a quasiclassical 
approximation, where the Fermi wavelength is much smaller 
than all other length scales. The Eilenberger equation reads: 32 

[Ef 3 + A, G] + iv F • VG - [E, G] = 0, (1) 

in which Vp is vector Fermi velocity of the quasi-particles in 
the superconducting region and S is a self-energy term related 
to e.g. elastic and spin-flip scattering centers. In the contrast 
of real part of £ which is an oscillatory function of energy, 
the imaginary part of the term is dissipative and decays when 
increasing the energy. Here, A and f 3 are defined as: 

^(l-t)- A =(-o)^=Gi)- 

H while □ and □ stand for 2 x 2 and 4x4 matrix quantities. 

In experimental situations, the diffusive regime of transport 
is often reached as very clean (ballistic) samples may be hard 




FIG. 1: (Color online) A): The schematic setup of a ferromagnetic 
trilayer of Ho|Co|Ho. The spiral curves show the trajectory of the 
magnetization vector in the Ho layer which rotates along a conical 
profile in the indirection, i) shows the configuration in which the 
magnetization patterns of the two Ho layers are completely identical, 
whereas in ii) the magnetization patterns follow a continuous spi- 
ral magnetization in the two Ho layers. B): Schematic model of the 
experimental setup of our proposed ferromagnetic Josephson junc- 
tion including a domain-wall and homogeneous ferromagnetic layer 
which enables a controllable triplet-supercurrent. The angle 6 rep- 
resents the orientation of the homogeneous magnetization in the Fl 
layer with respect to the z-direction. 

to fabricate. In the diffusive limit, impurities in specimen are 
very strong and consequently the self-energy term £ in the 
Eq. (2) dominates. The impurities will lead to randomization 
of the quasiparticle trajectories, scattering them all over k- 
space. By expanding the Green's function with respect to s- 
and p-wave spherical harmonics and performing an angular 
averaging process, one arrives at the Usadel-equation 28 : 

V(GVG) + -^{Ep 3 + dmg[h ■ <r,(h ■ <t) t },G] = 0, (2) 

where Eth — D/dp 2 is Thouless energy in which D is dif- 
fusive constant and dp is length of ferromagnetic layer, h is 
exchange field of ferromagnetic region and p 3 , a_ are Pauli 
matrixes which are available in Appendix A. 

Throughout the paper, we shall assume that the ferromag- 
netic layer has been sandwiched between two conventional 
s-wave superconducting leads whose interfaces are located at 
x = — dp/ 2 and dp/2. Due to the isotropic superconduct- 
ing order parameter and the impurity scattering, we may cap- 
ture all the essential physics by considering an effective one 
dimensional system, thus V = d/d x = d x in Eq. ( ). For in- 
vestigating the charge-current flowing through the system, we 
employ the following boundary conditions at the two contact 
regions with the superconducting reservoirs: 
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2(Gd x G= [GbcsWA 
2(Gd x G=[-G BCS (~cb),G} 



(3) 



where Gbcs is the bulk Green's function in the superconduc- 
tors and £ is defined as ratio between the resistance of the bar- 
rier region (Rb) and the resistance in the ferromagnetic film 
(Rf)- We disregard here the influence of spin-dependent in- 
terfacial phase-shifts occuring at the interfaces since their ef- 
fect is unimportant in the present context of an intrinsically in- 
homogeneous magnetization structure (including them would 
introduce slight shifts to the 0-ir transition points). 46 ' 47 . 

For solving the Usadel equation and implementing bound- 
ary conditions numerically, it is convenient to parameterize 
the Green's function. There are two standard parameteriza- 
tions approaches; 6- and Ricatti-parameterizations, and we 
will here employ the latter. The parameterized Green's func- 
tion then reads as follows: 



G = 



N(l - 77) 2N_1_ 
2N_j N(-l + 77) 



(4) 



By imposing a normalization condition for the Green's func- 
tion, namely G 2 — 1, N_ and N_ are obtained as 



7V = 



(5) 



1 + 77 I + 77 

Within the Ricatti-parametrization scheme, the components of 
the bulk superconductor Green's function are: 

7BCs(0) = ^/(l + c)e^/ 2 , 

7bcs0) = 7BCs(0)e" 10 , (6) 

where <fr is superconducting phase difference between the two 
s-wave superconducting leads and s, c is defined as sinh i? and 
cosh$, respectively, in which d = atanh(A(T)/£'). We use 
standard BCS temperature dependent of superconducting gap 
in our calculations and Ao = A(0) stands for superconduct- 
ing gap in the absolut zero. Throughout the paper we normal- 
ize all energies with respect to superconducting gap at the zero 
temperature (A ) and all lengths with respect to ferromagnetic 
layer length. We use units so that h = ks = 1. 

For investigating the electronic transport properties of all 
configurations one needs to obtain the Keldysh block of the 
Green's function. Under equilibrium conditions, the Keldysh 
component can be obtained from the Retarded and Advanced 
blocks as G K = (G R - G A )tanh(E/2k B T) and G A = 
— (p$G p$y . The charge-current is then obtained via: 



Ic_ 
loc 



dE Tr h 3 (dd x d) 



K 



(7) 



transport properties of several experimentally accessible con- 
figurations, in particular pertaining the experiment by Robin- 
son et a/. 35 . We underline once more that our approach grants 
us access to the full proximity regime in addition to compli- 
cated magnetization textures which cannot be treated analyti- 
cally. 



3. SPIN-TRIPLET SUPERCURRENT THROUGH A 
FERROMAGNETIC HO|CO HO TRILAYER 

In this section, we present main results of the paper: a the- 
oretical investigation of a ferromagnetic Ho|Co|Ho trilayer 
sandwiched between two s-wave superconductors, as recently 
experimentally studied in Ref. 35. The magnetization struc- 
ture of isolated Ho is experimentally known 43 and depicted in 
portion A) of Fig. . Analytically, the instantaneous direction 
of the local magnetization can thus be written as: 



h = h (cos a x + 5 sin a) 
{sin(Qx) y + cos(Qa;) z} , 



(8) 



where Q = 2tt/X and A is the spiral length in Ho. In Ref. 
35, this was estimated to A ~ 3.4 nm. The apex angle is 
denoted by a and equals 47r/9. 35 ' 43 For both Ho and Co fer- 
romagnetic layers, the strength of the exchange field is larger 
than alloys and compounds of Pd and Ni. For Pd-cNii-j; al- 
loys, weak strengths of the exchange field in the range of 
h/Ao sa 5 — 10 are accessible in the experiments. This should 
be contrasted with intrinsically ferromagnetic materials whose 
exchange field strengths are often very large (100 meV - 1 
eV). We set superconducting coherent length as £5 = 15 nm 
which is accessible for example in Nb. For more stability in 
the numerical code, we add an imaginary part to the quasi- 
particle energies equal to 8/Aq = 5 x 10~ 2 which model- 
ing inelastic scattering in the specimen. To simulate numeri- 
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in which Iqc = NobD/WcLf, Nq is the normal density of 
states per spin. Above, (^Gd x Gj denotes the Keldysh com- 
ponent of the Gd x G matrix. We now proceed to study the 



FIG. 2: (Color online) The normalized critical charge current 
for three types of magnetization textures. Top, middle and bot- 
tom frames show critical supercurrent through sandwiched uniform, 
domain-wall and spiral ferromagnetic layers, respectively. 
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FIG. 3: (Color online) The normalized critical supercurrent through 
a magnetic Ho|Co|Ho trilayer vs. the length of the Co layer for three 
values of Ho-layer lengths, Lho=2 nm, 5 nm and 8 nm. 
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FIG. 4: (Color online) The normalized critical supercurrent through 
a ferromagnetic trilayer of Ho | Co | Ho vs. the length of the Ho layer 
for three values of Co-layer lengths, Lco=2 nm, 5 nm and 8 nm. 



cally the ferromagnetic trilayer sandwiched between the two 
s-wave superconducting leads, we should first note the pa- 
rameter regime in which the quasi-classical Green's function 
method is valid. Due to the requirement of the Fermi energy 
being much larger than all energy scales, we set h/Ao = 80 
to model a strong exchange field still within the regime of va- 
lidity. Throughout our calculations, we fix the temperature at 
T/T c = 0.2 and also set the ratio of barrier and ferromagnetic 
resistances to ( — Rb / Rf = 5. For the value of the exchange 
field considered here, it can be shown that the spin-dependent 
phase-shifts occurring at the interface may be neglected. 

To begin with, we mention briefly the qualitative difference 
between having a homogeneous and inhomogeneous magne- 
tization in the magnetic layer. In Fig. , we show the critical 
charge current behavior vs. the thickness of ferromagnetic 
layer for three scenarios: a homogeneous exchange field, a 
domain-wall ferromagnet, and finally a spiral (conical) mag- 
netization texture. As seen, the fundamental difference be- 
tween the homogeneous case in the upper panel of Fig. 2 
and the two inhomogeneous scenarios is that the current be- 
comes long-ranged in the latter cases, i.e. the critical super- 
current decays on a much larger length-scale compared to the 
homogeneous case. This is due to the generation of a long- 
range triplet supercurrent which is sustained by the inhomo- 
geneous field 214 ' 23,40 ' 42 . Moreover, the spiral magnetization 
pattern gives rise to a rapid oscillation pattern superimposed 
on the 0-7T transitions, as was shown in Refs. 33,34. 

Now we model the ferromagnetic Ho | Co | Ho trilayer by: 

{h (cos ai + S sin a) x < — 
hz <x < ^ , (9) 

h (cos ai + S sin a) x > 

where we assume that the middle of the Co-layer is located at 
x = 0. In Fig. , we show the variations of the critical super- 



current when the Co-layer length is varying [see portion A) of 
Fig. for the structure under consideration]. Results are pro- 
vided for three distinct values of the Ho-layers length. We first 
assume that the Ho-layers have identical spiral magnetization 
patterns. As we shall see later, the results for the critical cur- 
rent are sensitive to the exact magnetization texture in the Ho- 
layers. The magnetization of the Co-layer is taken to be along 
the z-axis, i.e. parallel with the contact-interfaces, as should 
be reasonable for a thin-film structure. As seen in Fig. 3, for 
thin Ho-layers, the critical supercurrent displays only one 0- 
7T transition point. For larger values of Ln Q , the supercurrent 
decays monotonically similar to flowing critical supercurrent 
through an S|N|S junction where only spin-singlet condensa- 
tion contributes to the supercurrent. This is consistent with 
the experimental observation by Robinson et a/. 35 . For thick 
Ho-layers, the charge supercurrent thus behaves as if a normal 
layer has been sandwiched between two s— wave supercon- 
ducting. The long decay length of the supercurrent, compara- 
ble to a normal Josephson junction, is evidence of precisely a 
long-ranged spin-triplet supercurrent flowing through the fer- 
romagnetic Ho | Co | Ho trilayer. 

Next, we investigate how the current behaves upon alter- 
ing the Ho-layer thickness. We show results for three distinct 
values of Leo in Fig. 4. The Co-layer exchange field is as be- 
fore assumed to be oriented parallell to the interface regions. 
As seen, the current now decays in a non-monotonic fashion. 
Interestingly, Robinson et al? 5 observed a set of anomalous 
sharp peaks in the current when Lho is increasing. Although 
we are not able to reproduce such sharp peaks within this qua- 
siclassical treatment, we confirm the non-monotonic depen- 
dence of the critical current observed by Ref. 35. In order 
to investigate further if the specific magnetization profile is 
crucial with regard to the appearance of the anomalous behav- 
ior observed experimentally, we investigate a slightly different 
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FIG. 5: (Color online) The normalized critical supercurrent through 
a trilayer of Ho|Co|Ho vs. the length of the Ho layer for two values 
of Co-layer lengths, Lco=2 nm and 8 nm. In this case, we assume 
that the magnetization vector in Ho layers follow a continuous spiral 
pattern without being interrupted by the Co layer. 
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FIG. 6: (Color online) The normalized critical supercurrent through 
sandwiched Bloch-wall|Co|Bloch-wall trilayer vs length of Bloch- 
wall layers for three different values of Co layer lengths, Lc =2 nm, 
5 nm and 8 nm. Here magnetization direction of Co layer is along 
the z-axes. 



magnetization texture model in the Ho-layers. Whereas the 
magnetization pattern previously was assumed to be identical 
in both layers, we show in Fig. results for the case when 
the magnetization pattern in the right Ho layer couples con- 
tinuously to the left layer as if the Co-layer has no influence 
on it [see part z'z) of portion A) of Fig. ]. Unlike the first 
case, several minima now appear in the critical supercurrent, 
out of which two are 0-ir transition points while the others are 
local minima. This demonstrates that the exact behavior of 
the supercurrent is sensitive to the specific magnetization pro- 
file in the inhomogeneous magnetic layers. Motivated by this 
finding, we explore in the next section how the results change 
when the Ho-layers are replaced by domain-wall ferromag- 
nets. 



4. SPIN-TRIPLET SUPERCURRENT THROUGH A 
FERROMAGNETIC STRUCTURE WITH DOMAIN-WALLS 

We keep the same parameters as used in the previous sec- 
tion, but now replace the spiral magnetization patterns in por- 
tion A) of Fig. with domain-wall ferromagnets. Both a Neel 
and Bloch-wall configuration have been investigated by us nu- 
merically, and were found to give virtually identical results 
for the charge-current transport. Thus, we here present results 
only for the Bloch domain-wall case. Fig. 6 shows the crit- 
ical charge-current through a ferromagnetic trilayer with the 
following magnetization structure: 

!h(cos f3 y + sin fJ z) x<-^f 
hz <K^f (10) 

h(cos (3 y + sin j3 z) x > 



where dw is width of domain-wall and: 

/? = -2atanhf - T ~ Ld/2 y (11) 
V dw ' 

We set dw=L>d/2 throughout our computations, and provide 
results for the critical supercurrent biased through the trilayer 
vs. the domain-wall— layer length Ld for three distinct values 
of the Co-layer lengths in Fig. . For Lc = 2 nm the crit- 
ical supercurrent features two Q-tt transition points and these 
transition points disappear when Lq increases. This finding 
coincides with the Ho | Co | Ho trilayer. As can be seen from 
Fig. 6, the variation of critical supercurrent vs. the length 
of domain-wall— layer shows a non-monotonic behavior, as in 
Fig. 4. In effect, domain-wall ferromagnets can serve a simi- 
lar purpose as Ho with regard to the generation of a long-range 
current. 

Due to the somewhat complicated magnetization texture 
in the trilayer structure considered in Ref. 35, it is tempt- 
ing to consider if it is possible to simplify the structure 
of the magnetic layers and still obtain a setup where the 
triplet supercurrent can be experimentally controlled. In 
Ref. 15, a magnetic trilayer consisting of homogeneous, 
misaligned ferromagnets was proposed as a setup where the 
long-range current could be controlled by varying the an- 
gle of misalignment. Experimentally, it would nevertheless 
be highly challenging to exert individual control over the lo- 
cal magnetization field in each layer via application of exter- 
nal fields. We here propose another type of heterostructure 
which might be more beneficial in this regard. As is seen 
in part B) of Fig. , we consider a superconductor|domain- 
wall|ferromagnet|superconductor junction where the coercive 
field of the homogeneous ferromagnet is sufficiently low to 
allow tuning of its magnetization via an external field with- 
out altering the domain-wall ferromagnet. In this way, the 
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FIG. 7: (Color online) The normalized critical supercurrent through 
our proposed B loch-wall | homogeneous ferromagnetic bilayer vs. the 
length of homogeneous ferromagnetic layer Lpi for four values of 
magnetic orientation angle of homogeneous ferromagnet with re- 
spect to i-direction, 8=0, 7r/6, n/3 and 7r/2. 



domain-wall ferromagnet serves as a source for the long-range 
triplet current, while the orientation of the homogeneous fer- 
romagnet can tune this contribution. We define the magnetic 
field orientation angle of the Fl layer as 9 (see Fig. , portion 
B)), and thus the magnetization profile of the ferromagnetic 
layer reads: 



5. SUMMARY 

We have investigated the possibility of establishing a long- 
range spin-triplet supercurrent through an inhomogeneous fer- 
romagnetic region consisting of a Ho|Co|Ho trilayer sand- 
wiched between two conventional s-wave superconductors. 
Utilizing a full numerical solution in the diffusive regime of 
transport, the behavior of the supercurrent in several relevant 
configurations of the magnetic trilayer has been obtained. We 
find qualitatively very good agreement with the recently re- 
ported experimental results by Robinson et al. 35 regarding the 
behavior of the supercurrent as a function of the width of the 
Co-layer. Moreover, we find a synthesis of 0-ir oscillations 
with superimposed rapid oscillations when varying the width 
of the Ho-layer. We are not able to reproduce the anomalous 
peaks observed experimentally in this regime, but note that the 
results obtained are quite sensitive to the exact magnetization 
profile in the Ho-layers. We also investigate the supercurrent 
in a system where the intrinsically inhomogeneous Ho fer- 
romagnets are replaced with domain-wall ferromagnets, and 
find similar behavior as in the Ho|Co|Ho setup. In addition, 
we propose a novel type of ferromagnetic Josephson junc- 
tion involving a domain-wall and homogeneous ferromagnet 
which could be used to obtain a controllable spin-triplet su- 
percurrent. The advantageous of our proposed structure com- 
pared to the structure considered in Ref. 15 is a simpler mag- 
netic profile which could be beneficial from an experimental 
point of view. 



, _ f hi (cos 9 z + sin 9 y) 
1 /i2 (cos p y + sin j3 z) 



X < Lrp\ 

sin/3 z) x > L 



'Fl 



(12) 



As before, we set the exchange field in the homogeneous 
ferromagnetic region to h%/Ao = 15 whereas it is stronger 
in the domain- wall region hi/ 'Ao = 70 to avoid influence of 
the magnetic field on the domain-wall region. The width of 
the domain-wall is here set to — 3.4 nm, and we show the 
results for the critical supercurrent in Fig. 7 for four values 
of magnetic field orientation angle of the Fl layer, #=0, ir/6, 
7r/3, and tt/2. It it seen that the current is enhanced from 9=0 
up to values approximately near 9=%/6, suggesting that the 
triplet contribution to the current is tuned. The reason for this 
is that a singlet supercurrent would be completely invariant 
under a rotation of the exchange field since such a current is 
spinless. By increasing value of the angle further from 9=ir/6 
up to 7r/2, the critical supercurrent is suppressed. The main 
advantage of this setup compared to e.g. the trilayer structure 
considered in Ref. 15 is that only one ferromagnetic layer 
needs to have its magnetization orientation tuned, which is 
experimentally more feasible than individually controlling the 
magnetization structure of each individual layer in a trilayer 
structure. 
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Appendix A: Pauli matrixes 



The Pauli matrices we use in this paper 



1 



-i 



7-1 = 1 1 o >^ = i o '^=io -ir 



1 


1 = ' l) ' 1 = (o l)'^ = \i nl 



7~1 



o -in 



Pl =\n I r P2= Vn (rr p 3-\o -l 



1 
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